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Abstract 

Isospin violating valence and sea distributions are evaluated due to QED leading 
0(a) corrections to the standard QCD evolution equations. Unique perturbative 
predictions are obtained within the radiative parton model, and confirm earlier 
results. Nonperturbative contributions have been estimated and depend on a single 
free parameter chosen to be the current quark mass. The relevance of our predictions 
for extracting sin 2 9w from DIS v(p)N data ('NuTeV anomaly') is discussed as well. 



The NuTeV collaboration recently reported [1] a measurement of the Weinberg angle 
s^y = sin 2 9w which is approximately three standard deviations above those presented 
in [2] for the world average of other electroweak measurements. Possible sources for this 
discrepancy (see, for example, [3, 4, 5, 6, 7]) include, among other things, an isospin 
violating contribution 5RJ to the relation [8] Rp W — \ — s 2 ^, 
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This contribution is given, for N = \{j> + n), by [3, 5, 7] 
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where the second moments 5Q V (Q 2 ) and Q V (Q 2 ) of the valence distributions (5)q v (x, Q 2 ) 
are Q V (Q 2 ) = x q v (x,Q 2 )dx, q v = q — q, and 5Q V (Q 2 ) = x 5q v (x,Q 2 )dx with 

5u v (x,Q 2 ) = u p v (x,Q 2 )-d n v (x,Q 2 ) 

5d v (x,Q 2 ) = d p v (x,Q 2 ) -<(x,Q 2 ). (3) 

NLO QCD corrections do not significantly [4, 6, 7] contribute to (1). Nonperturbative 
calculations [9, 10, 11] of 8q v in (3) were found [5] to possibly reduce the discrepancy 
between the neutrino and electroweak measurements of s 2 ^ by up to 40%. A somewhat 
different calculation [12], based on QED contributions to 5q v (x, Q 2 ), resulted in a compa- 
rable reduction of this discrepancy. 

We shall approach the issue of the isospin violations within the framework of the 
radiative parton model [13] and obtain predictions for Sq v which depend on a single free 
parameter required by the nonperturbative contribution and chosen to be, as in [12], the 
current quark mass. Although our method differs from those in [5], the results turn out 
to be similar. On the other hand, our starting position is only a little different to [12], 
the results, as anticipated in that work, turn out to be quite similar. 



Following [12] we shall evaluate the modifications of the parton distributions due to 
QED radiative effects. In contrast to [12, 14] we shall calculate these effects only to 
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leading order in a, which is sufficient for our purpose here and furthermore simplifies 
the calculations considerably. The inclusion of QED 0{a) corrections modifies the QCD 
evolution equation for charged parton distributions by an additional term which, in an 
obvious symbolic notation, reads [15, 16] 

Q 2 ) = %;[P«*Q + P i9 *9] + ^ e 2 q PJ q * q (4) 

with P^ q (z) = (izf ) • Notice that the addition [12, 14] of a further term (a/2n) e 2 P 97 * 7 
to the r.h.s. of eq. (4) would actually amount to a subleading 0(a 2 ) contribution since 
the photon distribution 7(2;, Q 2 ) of the nucleon is of O(o) [17, 18, 19, 20, 21, 22]. The 
standard QCD evolution equation for the gluon distribution remains unaltered at leading 
O(a). The perturbative (p) 0(a) modifications of the valence quark distributions in (3), 
as implied by (4), are now given by 

5 p u v (x,Q 2 ) = I ding 2 f — P ( u v (y, q 2 ) 
2tt Jq2 J x y \yj 

5 P d v (x,Q 2 ) = [ Q d\nq 2 p(-)d v (y,q 2 ) (5) 

2tt J Q 2 J x y \yj 

with P(z) = (el - e^)P^ q (z) and Ql = /J[ = 0.26 GeV 2 , u v (x,q 2 ), d v (x,Q 2 ) are taken 

from the dynamical (radiative) parton model [13]. Notice that this treatment of the 

perturbative components differs from that in [12]. 

For the nonperturbative (np) modifications we estimate, following eq. (12) in [12], 

S np u v (x) = TT" ln% / — P (-] u v(Vi Ql) 
2vr i4 J x y \yj 

S np d v (x) = ~ ln^ f^i p(*) dv (y,Q 2 ) (6) 
2tt Hq J x y \y J 

where we take /x = rn q ~ 10 MeV, i.e. of the order of the current quark masses [2]. Here 
the parton distributions at Q 2 < Ql were taken to equal their values at the perturbative 
input scale Q 2 = Qjj, i.e. were 'frozen'. 

In fig. 1 we show the total 5 p q v (x, Q 2 ) + S np q v (x) and the purely perturbative 
S p q v (x,Q 2 ) isospin violating distributions at a typical scale of Q 2 = 10 GeV 2 . Our total 
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isospin violating distributions are quite similar to those in [12], as well as to those in [5] 
which were obtained by rather different methods. The fact that only our total results 
are compatible with nonperturbative (bag) model expectations [5, 9, 10, 11] is indicative 
for the necessity of our nonperturbative estimates in (6) and that the rather marginal 
perturbative contributions in (5) are not sufficient for a realistic estimate of isospin vio- 
lating effects. It is, furthermore, particularly interesting to note that our results depend, 
as mentioned above, only on a single free parameter, i.e. /x in (6). 

Encouraged by this agreement we present in fig. 2 and 3 the corresponding predic- 
tions for the isospin violating distribution 8u and 5d of sea quarks as defined in (3) and 
obtained from (5) and (6) by replacing u v and d v by u and d, respectively. At Q 2 = 10 
GeV 2 the perturbative contribution S p q(x, Q 2 ) in fig. 2 does not dominate the total result 
5 p q(x, Q 2 ) + S np q(x, Q 2 ) which is obviously in contrast to the predictions at higher scales 
as illustrated, for example, at Q 2 = in fig. 3. Similar results are obtained for the LO 
CTEQ4 parton distributions [23] where also valence-like sea distributions are employed 
at an input scale Ql = 0.49 GeV 2 , i.e., xq{x, Ql) — > as x — > 0. Such predictions may be 
tested by dedicated precision measurements of Drell-Yan and DIS processes employing 
neutron (deuteron) targets as well. 

As a possible immediate application of our predictions for the isospin violating valence 
distributions let us finally turn to the relation in (1). Since the isospin violation generated 
by the QED 0(a) correction is such as to remove more momentum from up-quarks than 
down-quarks, as is evident from fig. 1, it works in the right direction to reduce the NuTeV 
anomaly [1], i.e. sin 2 9 W = 0.2277 ± 0.0013 ± 0.0009 as compared to the world average 
of other measurements [2] sin 2 9w = 0.2228(4). This effect is also slightly Q 2 dependent 
because of the perturbative contribution in (5). At Q 2 = 10 GeV 2 , appropriate for 
the NuTeV experiment, our total distributions in fig. 1 imply 5U V = —0.002226 and 
SD V = 0.000890 which, together with U v + D v = 0.3648 in (2), leads to a change in the 
measured value of sin 2 9 W of S sin 2 9 W = SRJ = —0.0020 according to (1). Since the value 
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of sin 2 9w from the NuTeV experiment, with 5RJ = in (1), is 0.0049 larger than the 
world average of other measurements, our predicted isospin violation SRJ reduces this 
discrepancy by about 40%. This reduction may be slightly diminished if it is corrected 
for experimental acceptance cuts [3]. 

To summarize, we evaluated the modification of parton distributions due to QED lead- 
ing 0(a) corrections to the standard QCD evolution equations. For the isospin violating 
valence Sq v (x,Q 2 ) and sea 5q(x,Q 2 ) distributions (q = u, d) unique perturbative predic- 
tions are obtained within the dynamical (radiative) parton model. The nonperturbative 
contributions to 5q v and 5q have been estimated and depend on a single free parameter 
chosen to be the current quark mass. The results for the valence distributions bear some 
similarity to nonperturbative bag-model expectations. Our total predictions for 5q v re- 
duce significantly the discrepancy between the large result for sin 2 6 W as derived from deep 
inelastic v{v)N data ('NuTeV anomaly') and the world average of other measurements. 
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Figure 1: The isospin violating 'majority' 5u v and 'minority' Sd v valence quark distribu- 
tions at Q 2 = 10 GeV 2 as defined in (3). The perturbative GRV(pert.) predictions are 
calculated according to (5). With the non-perturbative contribution in (6) added, one 
obtains the total estimates GRV(total). The bag model estimates are taken from ref. [5]. 
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Figure 2: The isospin violating sea distributions 8u and 5d at Q 2 = 10 GeV 2 as defined 
in (3) with u v , d v replaced by u, d. The perturbative (pert.) predictions are calculated 
according to (5) with u v , d v replaced by u, d. The same replacement holds for (6) when 
calculating the non-perturbative contributions which have to be added to the 'pert.' 
predictions in order to obtain the total results. The LO GRV sea distributions [13] are 
used throughout. 
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